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We  report  the  observation  of  resolved  atomic  interaction  sidebands  (ISB)  in  the  87  Sr  optical  clock 
transition  when  atoms  at  microkelvin  temperatures  are  confined  in  a  two-dimensional  optical  lattice.  The 
ISB  are  a  manifestation  of  the  strong  interactions  that  occur  between  atoms  confined  in  a  quasi-one- 
dimensional  geometry  and  disappear  when  the  confinement  is  relaxed  along  one  dimension.  The 
emergence  of  ISB  is  linked  to  the  recently  observed  suppression  of  collisional  frequency  shifts.  At  the 
current  temperatures,  the  ISB  can  be  resolved  but  are  broad.  At  lower  temperatures,  ISB  are  predicted  to 
be  substantially  narrower  and  useful  spectroscopic  tools  in  strongly  interacting  alkaline-earth  gases. 
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Experimental  efforts  in  control  of  ultracold  alkaline- 
earth  atoms  (AEA)  have  led  to  remarkable  developments 
in  optical  atomic  clocks  [1,2],  which  are  approaching  the 
accuracy  of  single  ion  standards  [3].  Fermionic  AEA  are 
also  attracting  considerable  attention  in  the  context  of 
quantum  information  science  [4-6].  Many  of  these  appli¬ 
cations  require  reaching  the  strongly  interacting  regime  in 
which  interactions  dominate  over  all  relevant  dynamical 
energy  scales.  In  this  Letter,  we  demonstrate  the  onset  of 
strong  interactions  in  a  low  dimensional  gas  of 87  Sr  atoms 
by  observing  resolved  interaction  sidebands  (ISB)  in  Rabi 
spectroscopy  of  the  lS0-3P0  optical  clock  transition. 

Rabi  spectroscopy  has  served  as  an  excellent  probe  for 
interacting  alkali  gases.  For  example,  it  was  used  to  char¬ 
acterize  Bose  Einstein  condensation  (BEC)  in  spin- 
polarized  hydrogen  [7],  to  measure  fermion  pairing  in  the 
BEC-BCS  crossover  [8],  and  to  resolve  the  Mott  insulator 
shell  structure  in  bosonic  lattice  systems  [9].  Here,  we 
show  that  the  ISB  that  emerge  during  optical  Rabi  inter¬ 
rogation  can  become  a  powerful  spectroscopy  tool  to 
uniquely  probe  lattice-trapped  AEA.  At  microkelvin  tem¬ 
peratures  and  loading  conditions  of  mainly  one  or  two 
atoms  per  lattice  site,  we  observe  a  single  ISB,  which 
can  be  resolved  from  the  carrier  but  is  broad.  However, 
in  the  quantum  degenerate  regime,  the  ISB  are  expected  to 
become  substantially  narrower  and  could  thus  lead  to  a 
precise  determination  of  the  s-wave  15'0-3P0  interaction 
parameters  that  remain  unknown.  We  also  show  that  the 
variation  in  ISB  spectra  between  sites  with  different  num¬ 
bers  of  atoms  allows  for  a  precise  characterization  of  atom 
number  occupations  in  interacting  AEA  systems. 

When  the  lS0  and  3P0  electronic  degrees  of  freedom  are 
represented  as  an  effective  pseudospin  1/2  (g  and  e ),  two 
nuclear-spin-polarized  atoms  form  a  simple  model  for 
understanding  the  appearance  of  ISB  (Fig.  1).  We  consider 
atoms  trapped  in  a  tube  geometry  with  a  weak  trapping 
frequency  coz  along  Z  and  a  strong  transverse  confinement 
hcoXY  ^  kBTXY,  large  enough  to  forbid  any  transverse 


dynamics  [TX(TY)  is  the  temperature  along  X(F)].  We 
assume  the  atoms  are  prepared  in  the  pseudospin  state 
|gg)  and  interrogated  by  a  linearly  polarized  laser  with 
bare  Rabi  frequency  ClB  and  detuning  from  the  atomic 
resonance  <5.  The  Pauli  exclusion  principle  forces  atoms 
in  identical  internal  states  to  occupy  different  vibrational 
levels  along  Z,  n  =  (nh  n2),  since  their  spatial  wave 
function  must  be  antisymmetric.  Being  identical  fermions, 
the  atoms  initially  experience  zero  s-wave  interaction. 
If  the  atoms  are  coherently  driven  to  e  [i.e.,  if  they  expe¬ 
rience  the  same  Rabi  frequency  £1%  =  (flWi  +  12^/2],  the 
effective  spin  remains  symmetric  under  exchange  and 
during  the  excitation  the  atoms  will  not  experience  any 
s-wave  interactions.  We  refer  to  the  spin- symmetric  states, 
|gg)=  \S  =  \,M  =  -1),  \ee)=  |1, 1),  and  {\eg)+\ge))/ 
72  =  11,0),  as  the  triplet  states  [10,11]. 

However,  in  optical  transitions,  a  small  component  of 
the  probe  beam  along  Z  leads  to  a  slightly  different  Rabi 
frequency  kln.(r]z)  for  each  mode,  with  rjz  =  kza^0/\/2 
the  Lamb-Dicke  parameter,  m  the  atom  mass,  kz  the 


FIG.  1  (color  online).  Schematic  diagram  of  Rabi  spectros¬ 
copy  in  the  strongly  interacting  regime.  When  <5  —  0,  only  the 
noninteracting  triplet  states  are  accessible,  giving  rise  to  an 
interaction-free  carrier.  However,  when  8  ~  U^,  the  singlet 
becomes  resonant,  and  a  slightly  inhomogeneous  coupling 
(A  ^  0)  can  induce  population  transfer  to  it,  as  manifested 
in  an  ISB. 
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component  of  the  probe  laser  wave  vector  along  Z,  and 
a^Q  =  yjh/m(x)z  the  Z  harmonic  oscillator  length.  If 
An,  =  (nBl  -  n„2)/V2  is  not  zero,  then  the  optical 
excitation-induced  inhomogeneity  can  transfer  some  of 
the  atoms  to  the  interacting  singlet  state,  (|  eg)  — 
\ge))/'Jl  =  |0,  0),  separated  from  |1,  0)  by  an  interaction 
energy  Un. 

In  the  limit  where  S  —  a  small  cannot 

overcome  the  energy  cost  required  to  drive  the  transition  to 
the  singlet  (Fig.  1).  Evolution  into  the  singlet  is  blocked 
and  the  atoms  remain  noninteracting.  However,  if  we 
increase  the  detuning  to  8  —  Un,  the  transition  |1,  —  1)  — ► 
|0,  0)  becomes  resonant.  Here,  a  small  can  efficiently 
transfer  the  atoms  to  the  singlet,  giving  rise  to  an  ISB. 

The  ISB  resemble  motional  sidebands  that  develop  when 
the  trapping  frequency  is  much  larger  than  the  recoil 
energy,  i.e.,  the  Lamb-Dicke  regime  [12]  where  the  carrier 
becomes  free  of  Doppler  shifts.  In  our  case,  the  ISB  are 
pushed  away  from  the  carrier  due  to  the  large  energy 
separation  between  interacting  and  noninteracting  states, 
leaving  the  carrier  free  from  interaction  effects,  up  to  a 
correction  proportional  to  XFll/U^.  Operating  our  Sr  lat¬ 
tice  clock  in  the  strongly  interacting  regime  has  allowed  us 
to  suppress  the  collisional  frequency  shift  to  the  10“ 17 
level,  as  reported  in  Ref.  [1].  The  observation  of  ISB 
confirms  the  connection  between  suppression  of  collisions 
and  strong  interactions,  supporting  the  prospect  of  operat¬ 
ing  neutral  atom  clocks  with  the  highest  possible  stability 
and  accuracy. 

Our  experiment  employs  nuclear  spin-polarized  (m7  = 
+  9/2)  87 Sr  atoms  loaded  in  two  lattice  configurations:  a 
one-dimensional  (ID)  vertical  lattice  which  creates  an 
array  of  quasi-two-dimensional  traps  (pancakes)  along  Y 
with  weak  confinement  along  the  X  and  Z  directions,  and  a 
deep  2D  lattice  which  creates  an  array  of  quasi- ID  traps 
(tubes)  in  the  X-Y  plane  and  has  weak  confinement  along 
Z.  To  prepare  the  atomic  system,  we  laser  cool 87  Sr  atoms 
to  about  2  julK  inside  a  magneto-optic  trap  (MOT)  based  on 
the  weak  lS0-3P1  transition  and  then  load  them  into  a  ID 
vertical  lattice.  The  spatial  distribution  of  occupied  ID 
lattice  sites  is  determined  by  the  vertical  extent  of  the 
MOT  cloud,  which  is  approximately  Gaussian  with  a  stan¬ 
dard  deviation  av  =  30  jmm.  For  the  2D  lattice  we  then 
adiabatically  ramp  up  the  horizontal  lattice  (along  X).  To 
remove  any  atoms  trapped  in  the  ID  vertical  lattice  outside 
the  2D  intersection  region,  we  ramp  the  vertical  lattice  off 
and  then  back  on.  The  number  of  horizontal  lattice  sites 
occupied  is  determined  by  the  radial  temperature  of  the 
vertical  lattice.  Atoms  load  into  —100  rows  of  tubes  uni¬ 
formly  distributed  along  F,  while  the  columns  distributed 
along  X  are  loaded  according  to  a  Gaussian  distribution 
with  standard  deviation  crH  of  6-10  jmm .  After  forming  the 
lattices,  we  perform  Doppler  and  sideband  cooling  using 
the  lS0-3P1,  F  =  11/2  transition.  Simultaneously,  atoms 
are  optically  pumped  to  the  m7  =  +9/2  ground  state 


magnetic  sublevel,  using  cr+  -polarized  light  on  the 
1S0-3P1,  F  =  9/2  transition,  propagating  along  a  bias 
magnetic  field  parallel  to  Z.  We  perform  spectroscopy  of 
the  1S0-3P0  transition  using  a  narrow  linewidth  laser  prop¬ 
agating  along  Y. 

We  quantify  the  number  of  atoms  by  detecting  fluores¬ 
cence  on  the  strong  1S0-lP1  transition  at  461  nm. 
Temperature  and  trap  frequencies  are  extracted  from 
Doppler  spectroscopy  along  Z  and  vibrational  sideband 
spectra  [1,13].  We  estimate  that  we  loaded  —7000  atoms 
in  the  ID  lattice  distributed  in  —420  pancakes  and  —17 
atoms  per  pancake.  We  determine  coY  =  27t  X  80  kHz 
and  o)x  z  =  2i t  X  500  Hz,  and  TXYZ  ~  4  A&K.  For  the 
2D  lattice,  we  loaded  —3000  atoms,  and  —30%  of  the 
populated  tubes  are  doubly  occupied.  We  determine 
o)x  =  2tt  X  110  kHz,  o)Y  =  2tt  X  70  kHz,  and  (x)z  = 
2i t  X  800  Hz  at  trap  center,  and  Tx,y,z  ~  4-5  /ulK. 

ISB  develop  in  the  line  shape  only  in  the  parameter 
regime  where  y  =  Eint/(Cl)T  »  1.  Here,  Eint  =  (N  —  1)X 
\{U)t\/2  is  the  mean  interaction  energy  per  particle 
and  N  represents  the  number  of  atoms  per  lattice  site. 
(U)T  characterize  the  strength  of  two-body  interactions 
and  depends  both  on  the  confinement  volume  and 
temperature  as  [14]  (U)T  —  w#(^^)#(^;)#(^^)  with 

u  =  4 a~g^m(oxcoYcoz/ h ,  a~g  the  singlet  1S0-3P0  scatter¬ 
ing  length  [6]  and  #  and  #  functions  characterizing 
the  temperature  dependence,  #(v  »  1)  —  2#(v  »  1)  — ►  1 
and  #(v  «  1)  —  #(v  «  1)  — ►  y/x. 

Using  the  experimental  conditions  and  assuming  a  mod¬ 
erate  a~g  =  —  70a0  (a0  the  Bohr  radius),  we  obtain  y2D  — 
10  in  the  doubly  occupied  sites  at  the  center  of  the  2D 
lattice.  On  the  contrary,  for  the  central  pancake  of  the  ID 
lattice  with  a  mean  number  of  17  atoms  and  the  same 
scattering  length,  y  is  reduced  to  y1D  —  0.6.  Based  on 
these  values,  one  predicts  the  development  of  ISB  only  in 
the  2D  lattice  geometry,  a  prediction  confirmed  by  our 
measurements.  We  measured  a  series  of  line  shapes  in 
the  ID  and  2D  lattices.  In  the  2D  lattice  case  we  also 
varied  the  pulse  area  and  Rabi  frequency.  The  resulting 
data  are  shown  in  Fig.  2.  To  remove  some  of  the  statistical 
noise  we  performed  multiple  scans  across  the  transition 
resonance  under  each  experimental  condition.  For  each 
scan,  the  drift  of  our  ultrastable  laser  was  canceled  to 
<50  mHz  per  second,  and  the  laser  was  stepped  by  2  Hz 
each  experimental  cycle  (1.5  s).  The  center  of  each  scan 
was  determined  by  a  fit  to  a  Lorentzian.  Each  data  set 
consists  of  20  superposed  scans  and  the  direction  of  the 
scan  was  alternated  to  reduce  systematic  effects  due  to 
residual  drift.  We  separated  the  data  into  bins  and  calcu¬ 
lated  the  mean  and  standard  error  of  the  mean  for  each  bin. 
At  our  4-5  /ulK  temperatures,  the  ISB  in  the  2D  lattice  are 
noticeable  but  broad.  To  improve  resolution,  we  reflect  the 
binned  line  shape  at  positive  S  (free  of  ISB  when  a~g  <  0) 
and  subtract  it  from  the  corresponding  bin  at  negative  8 
(see  Fig.  4). 
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ISB  in  the  tightly  confined  geometry  can  be  understood 
quantitatively  by  considering  an  isolated  tube  with  N 
atoms.  We  characterize  the  tube  dynamics  by  defining  a 
set  of  effective  spin  operators,  S%J’Z,  in  the  {e,  g}  basis 
[1,11].  Here  the  subscript  rij  is  drawn  from  a  fixed  set  n  = 
{ni,  n2, . . nN}  of  initially  populated  vibrational  modes 
along  Z.  The  description  of  the  system  in  terms  of  effective 
spin  operators  is  valid  provided  those  initially  populated 
modes  remain  singly  occupied  during  the  excitation  pro¬ 
cess.  In  the  rotating  frame,  the  Hamiltonian  of  the  system 
up  to  constant  terms  becomes  [11] 


TjS 

"a  =  -ssz 


N 

X  aH]s*nj 

j= i 


£  Unj.n, 

f*j 


(1) 


Sz,x  =  X'/l  |  SnX  are  collective  spin  operators.  The  quantity 
U„.,n/  =  measures  the  strength  of 

the  interactions  between  two  atoms  in  the  singlet  state. 
Inpn,  is  an  overlap  integral  between  harmonic  oscillator 

modes  along  Z  [14].  Only  the  Z-mode  distribution  is 
treated  exactly  since  at  current  temperatures  only  a  few 
transverse  excited  modes  are  populated.  The  population  of 
transverse  modes  is  accounted  for  as  a  renormalization  of 
the  interaction  parameter. 

The  interaction  part  of  the  Hamiltonian  is  diagonal  in 
the  collective  spin  basis  \S,  M),  S  =  00, . . N/2  and 
\M\  <  S.  For  N  =  2,  the  spin  basis  is  spanned  by  the  triplet 
states  |1,  M)  and  the  singlet  |0,  0).  Among  the  collective 
states  only  the  S  =  N/2  states  are  noninteracting.  States 
with  S  <  N/2  experience  a  finite  interaction  energy.  In  the 
presence  of  excitation  inhomogeneity,  S  is  no  longer 


FIG.  2  (color  online).  Measured  Rabi  line  shapes  for  g  to  e 
interrogation  (normalized  with  respect  to  peak  height).  The  main 
panel  shows  line  shapes  taken  in  a  ID  lattice  (black  triangles) 
and  in  a  2D  lattice  (connected  green  circles)  when  the  transition 
is  excited  with  a  i t  pulse  with  80  ms  duration.  The  noninteract¬ 
ing,  homogeneous  Rabi  line  shape  for  an  80  ms  tt  pulse  is  also 
shown  (solid  orange  line).  The  inset  shows  line  shapes  measured 
in  a  2D  lattice  using  a  160  ms  2tt  pulse  (connected  red  dia¬ 
monds)  and  a  160  ms  tt  pulse  (connected  blue  squares).  Each 
line  shape  was  measured  under  the  trapping  and  temperature 
conditions  described  in  the  main  text.  The  data  from  each  line 
shape  are  collected  into  bins  1-5  Hz  wide. 


conserved,  and  during  excitation  of  the  clock  transition, 
atoms  can  be  transferred  mainly  between  S  =  N/2 
and  S  =  N/2  —  1  states.  The  N  —  1  collective  excitation 
modes  \S  =  N/2-  l,M  =  N/2  -  1,  q)  (q  =  1,...,N  —  1) 
with  energies  UiN  can  be  accessed  from  the  initially 
g-polarized  state  \N / 2,  ~N/2)  using  <5  ~  UfN.  They 
give  rise  to  N  —  1  sidebands  in  the  line  shape  (Fig.  3): 

m  «  Nf(t,  s,  A*)  +  ^  f(t,  s  -  uf,  Anf  *),  (2) 

q=\ 

with  fit,  8,  y)  =  pr^sin2(V>’2  +  <52/2)  and  AUf*  = 
2IAA/2*  N/2\Sxj\N/2  -  l,  N/2  -  1,  q).  The  first 
term  is  associated  with  the  carrier  and  the  second  with 
the  ISB.  For  the  case  N  =  2,  we  have  U\’2  =  U and 
A  Ill’2  =  At  quantum  degeneracy,  only  the  lowest 

lying  vibrational  modes  need  to  be  considered,  n  = 
{0,  1, . . N  —  1}.  In  this  case,  the  interaction  sidebands 
are  very  narrow.  Assuming  t  =  stt/ ClB,  with  s  a  numerical 
constant  that  describes  the  pulse  duration,  the  peak  width  is 
just  2 ClB/s.  Measuring  the  interaction  sidebands  thus  pre¬ 
cisely  determines  a~g  and  characterizes  the  S  =  N/2  —  1 
spectrum,  which  has  structure  for  N  >  2.  Since  the  peak 
height  of  the  interaction  sidebands  scales  as  rj |,  large 
probe  inhomogeneities  are  required  to  observe  them  at 
low  temperature. 

In  Ref.  [9],  interaction-dependent  transition  frequency 
shifts  in  rf  spectroscopy  were  used  to  distinguish  sites 
with  different  occupation  numbers,  revealing  the  shell 
structure  of  the  Mott  insulator  phase.  The  g  and  e  states 
in  AEA  have  no  hyperfine  structure,  and  hence  rf  spectros¬ 
copy  is  not  applicable.  However,  as  shown  in  Fig.  3,  the 
optical  analog  of  rf  spectroscopy — the  ISB  spectroscopy — 
can  distinguish  sites  with  different  atom  numbers.  This 
capability  can  have  important  implications  in  quantum 
information  science  [4,6]. 

At  finite  T,  instead  of  a  fixed  set  of  populated  modes,  a 
thermal  average  should  be  taken.  To  correspond  to  our 
current  experiment,  we  focus  on  the  N  =  2  case.  After 
evaluating  the  thermal  average  [14],  we  obtain  the 


6  (Hz) 

FIG.  3  (color  online).  Excitation  fraction  vs  8  at  T  =  0  for 
lattice  sites  with  N  =  2-5  atoms.  The  ISB  reveal  the  S  =  N/2  —  1 
spectrum  and  can  be  used  to  probe  occupation  number 
in  a  lattice.  The  parameters  used  for  the  plot  were  r]z  =  0.4, 
n5  =  2tt  X  5  Hz,  t  =  1.57 r/D5,  and  u  =  2tt  X  2800  Hz. 
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FIG.  4  (color  online).  Finite  temperature  ISB  in  a  2D  lattice 
( g  to  e ).  Symbols  correspond  to  experimental  data  obtained  after 
reflecting  the  binned  line  shape  at  positive  detuning  (free  of  ISB) 
and  subtracting  it  from  the  corresponding  bin  at  negative  detun¬ 
ing.  Data  conditions  previously  shown  in  Fig.  2  are  displayed 
with  the  same  symbols  and  colors.  The  data  indicated  by  purple 
crosses  correspond  to  an  80  ms  tt  pulse  but  with  different 
trapping  conditions  (cox,y,z  =  2tt  X  {92,  60,  0.6}  kHz  at  trap 
center,  Tx,y,z  ~  {2.5,  2.5,  7.2}  /xK).  The  inset  shows  the  absence 
of  ISB  in  a  ID  lattice.  Solid  curves  correspond  to  theoretical 
predictions.  The  dip  at  8  ~  —70  Hz  in  the  measured  ISB  is  due 
to  residual  rrij  =  7/2  atoms. 


following  approximate  expression  for  the  ISB,  valid  under 
the  conditions  |<5|  >>  flB,  kBTz  >  hcoz,  and  t  =  stt/£Ib  « 
which  are  satisfied  in  our  experiment: 


(Ne)T 


rj^v^sCl^  (kBTzy/(U)Ty 
2\(U)t\  \hcoz)  6 


(3) 


Figure  4  shows  comparisons  between  the  measured  ISB 
and  theoretical  curves  calculated  using  Eq.  (3).  In  the  2D 
lattice,  sites  near  the  wings  of  the  lattice  beams’  Gaussian 
intensity  profiles  are  significantly  shallower  than  near  the 
center  of  the  beam  intersection  region,  i.e.,  co j_(X,  Y)  ~ 
£y_i_(0,  0)  exp[— (X2  +  Y2)/Wx],  with  co  x  =  ^JcoxcoY  and 
W x  =  -s/WxWy  ~  30  jmm  the  mean  beam  waist.  This  im¬ 
plies  u  varies  from  tube  to  tube;  so  to  compare  with 
experiment,  it  is  important  to  average  Eq.  (3)  over  the 
populated  tube  distribution.  We  assume  a  uniform  distri¬ 
bution  of  doubly  occupied  tubes.  With  this  procedure,  we 
obtain  line  shapes  that  reproduce  the  experimental  data  if 
we  use  a~g  =  —  280a0  and  77 z  ~  0.07(0.05)  for  the  cox  = 
2i t  X  87(73)  kHz  cases.  The  value  of  a~g  that  fits  the 
experimental  data  is  outside  the  range  of  effective  scatter¬ 
ing  lengths  a~g - (35-50)a0  reported  in  Ref.  [1]. 

However,  the  discrepancy  is  attributed  to  the  fact  that 
a~g  in  Ref.  [1]  corresponded  to  the  effective  scattering 
length  required  to  model  the  entire  lattice  array  by  a  single 
tube  at  the  trap  center,  and  thus  \a~* |  should  be  smaller 
than  \a~g\.  In  agreement  with  theory,  the  ISB’s  position  and 
width  are  on  the  order  of  (U)T  and  their  resonant  peak 
heights  (normalized  to  the  carrier)  decrease  with  decreas¬ 
ing  flB  and  increase  with  pulse  area  sir. 

We  determine  the  nuclear-spin  purity  of  the  atomic 
sample  to  be  >97%  by  scanning  the  probe  laser  over  the 


clock  transition  frequencies  for  other  nuclear-spin  states. 
The  signature  of  the  remaining  3%  is  the  small  resonant 
peak  at  8  ~  70  Hz  in  Fig.  2.  By  comparing  the  amplitude 
of  this  small  peak  with  the  ISB,  we  conclude  that  the 
observed  sidebands  are  not  caused  by  v-wave  collisions 
between  unpolarized  atoms. 

In  order  to  investigate  the  effect  of  tunneling,  we  also 
took  one  set  of  data  under  conditions  which  reproduced 
those  used  in  Ref.  [1].  The  presence  of  ISB  in  this  line 
shape,  measured  with  TXY  ~  2.5  julK  (see  Fig.  4),  com¬ 
pletely  rules  out  the  possibility  that  these  spectroscopic 
features  are  caused  by  tunneling  since,  under  these  con¬ 
ditions,  we  predict  that  tunneling  rates  are  suppressed  by 
2  orders  of  magnitude  compared  to  TX  Y  ~  4.5  julK. 

The  studies  presented  here  demonstrate  that  tools  devel¬ 
oped  for  clock  experiments  can  provide  precise  under¬ 
standing  of  strongly  interacting  many-body  systems. 
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